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We examine the impact of the combination of a static electric field and a non resonant linearly 
polarized laser field on an asymmetric top molecule. Within the rigid rotor approximation, we an- 
alyze the symmetries of the Hamiltonian for all possible field configurations. For each irreducible 
representation, the Schrodinger equation is solved by a basis set expansion in terms of a linear 
. combination of Wigner functions respecting the corresponding symmetries, which allows us to dis- 

^SJ ' tinguish avoided crossings from genuine ones. Using the fluorobenzene and pyridazine molecules as 

prototypes, the rotational spectra and properties are analyzed for experimentally accessible static 
field strengths and laser intensities. Results for energy shifts, orientation, alignment and hybridiza- 
tion of the angular motion are presented as the field parameters are varied. We demonstrate that a 
proper selection of the fields gives rise to a constrained rotational motion in the three Euler angles, 
the wave function being oriented along the electrostatic field direction, and aligned in the other two 
angles. 
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I. INTRODUCTION 



The manipulation of large molecules by using external fields represents, in spite of its long history, a very active 
and promising research area. Indeed, major efforts have been undertaken to create samples of oriented and/or 
aligned molecules, and a large variety of experimental techniques have been developed, such as, e.g., the brute force 
orientation [l| , hexajjole focusing , a train of laser pulses d, Q , or a combination of a laser pulse and a weak 
static electric field 0-0- The control and manipulation of the directional features of molecules, i.e., of their rotational 
degree of freedom, optimize the information content on experimental measurements performed in the laboratory frame. 
f~| , Indeed, the availability of asymmetric top molecules in oriented and /or aligned pendular states allows for a wealth 
Qh' of interesting applications in areas as diverse as spectroscopy [l^, [HI , photoelectron angular distributions [l2l. [l3j , 
. sterodynamic control of chemical reactions [l3 - [T7| , dissociation of molecules [l8l - l2]| , electron diffraction |22l | , or 
' high-harmonic generation [2^, [23| • 
^ , The experimental achievements have been accompanied by theoretical efforts to understand and explain the in- 
■ triguing physical phenomena appearing in asymmetric top molecules exposed to external fields. Regarding the impact 
OO I of radiative fields on these molecules, the corresponding theoretical studies have been especially fruitful in explaining 
I a vast amount of experimental results, such as, e.g., the rotational revival structure following the irradiation by an 
' I ' intense picosecond laser pulses [l^ , the three-dimensional alignment by elliptically polarized laser fields j26l - l28| or the 
\G , use of long and short laser pulses to control the rotation [2^ [s^l ■ Analogously, the motivation of the theoretical works 
' considering an electrostatic field was either to interpret some experimental results fiol. o r to confirm the feasibility of 
' other experiments, e.g., the Stark deceleration of polyatomic asymmetric molecules |3l| . The molecular orientation 
due to the interaction with a static electric field has been investigated for asymmetric top molecules with their per- 
manent dipole moment /x parallel to a principal axis of inertia, and for the non-parallel case .32, . ,33j . In the strong 
electrostatic field regime, an analytical study of the energy-level representation has shown that the asymmetric top 
pendular states are well described by a two-dimensional anisotropic harmonic oscillator |34|. and it has been used to 
c5 , reproduce spectroscopic results in the pendular regime for static fields up to 200 kVcm~^ 

A detailed analysis of the rotational spectrum of symmetric top molecules exposed to combined electrostatic and 
nonresonant radiative fields was recently performed by Hartelt and Friedrich jssj . For tilted fields, only the projection 
of the total angular momentum J onto the body fixed frame z-axis K remains as good quantum number, and a 2D 
description of the rotational spectrum of the molecule is required. The corresponding dynamics is very complicated, 
indeed, in the presence of a static electric field, it has been shown that the molecular spectrum presents classical and 
quantum monodromy 13 6| . They provide correlation diagrams between the field- free states and the pendular levels of 
the intense laser field [33 as well as the strong electrostatic field regime. For a selection of states, they investigate 
the energy shifts and directional properties (orientation and alignment) for parallel and perpendicular fields. In these 
systems, the coupling of both field interactions could provoke an enhancement of the orientation giving rise to an 
oriented and antioriented pair of levels. For an oblate system, this phenomenon appears in the tunneling doublets 
created by the interaction of the molecular polarizability with the linear polarized laser field (this effect was already 
analyzed for linear molecules [38l - l40j ). whereas for a prolate molecule, it appears among exactly degenerate doublets 
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of indefinite parity appearing in the strong laser field regime. 

A classical theoretical analysis of asymmetric top molecules exposed to a combination of static and laser fields has 
been performed recently [4l|. However, the quantum analog has not yet been addressed in the literature, to the best 
of our knowledge. Recently, the authors have developed a diabatic model to describe the evolution of alignment and 
orientation of asymmetric top molecules in combined fields as the laser intensity is varied (4^ . The outcome of this 
theoretical study has been compared to the experimental data obtained for the benzonitrile molecule [l3l | proving 
the importance of non-adiabatic processes in the field-dressed molecular dynamics. Thus, motivated by the current 
experimental interest on these asymmetric molecules [3iEj 12, 13j and by the fact that the rotational dynamics of 
most polyatomic molecules can be described as asymmetric tops, we extend in the present work the previous study 
on symmetric tops [ssf to these more complicated systems. We perform a theoretical investigation of an asymmetric 
top in the presence of combined electrostatic and nonresonant radiative fields within the rigid rotor description. 
The field-dressed rotational spectrum is significantly more complicated, and the more general case of no collinear 
field requires a full 3D description. We will perform a detailed analysis of the symmetries of the Hamiltonian for 
all possible field configurations. In tilted fields, the reduction of the symmetries enhances the complexity of the 
spectrum, and a large amount of avoided crossing appears between states of the same symmetry. Hence, to simplify 
the analysis and interpretation of our results, the Schrodinger equation is numerically solved for each irreducible 
representation by expanding the wave function in a basis with the corresponding symmetry. As prototype examples 
we consider (CgHsF) and pyridazine (C4H4N2) molecules. These two systems have similar values of their polarizability 
tensors and dipole moments, but different inertia tensors and are, therefore, affected differently by the external fields. 
We explore their rotational spectrum as either the laser intensity, the electrostatic field strength, or the inclination 
angle between them is varied. Our focus is on the energy shifts, the directional properties and the hybridization of 
the angular motion. Depending on the dominant interaction, a rich field-dressed dynamics is observed with levels 
achieving different degrees of orientation and/or alignment. The role played by the inclination angle is exemplary 
investigated via a set of states and in avoided crossings between two adjacent levels. Moreover, we show that due to 
the combination of both field interactions the rotational motion is restricted in the three Euler angles, being oriented 
along the static electric field direction and constrained in the XY plane of the laboratory frame, which is perpendicular 
to the laser polarization. This mechanism of orientation and 2D alignment is very sensitive to the field parameters 
and to the molecular properties. 

The paper is organized as follows. In Sec. |lT]the rotational Hamiltonian is presented together with a comprehensive 
consideration of its symmetries for the different field configurations. In Sec. IIIIl wc discuss the numerical results 
for two asymmetric molecules, fluorobenzene and pyridazine, as the field parameters arc modified. In particular, we 
explore three different cases: i) for fixed laser intensity and three inclination angles, we vary the electrostatic field 
strength; ii) for fixed electrostatic field and three inclination angles the laser intensity is enhanced; and iii) for fixed 
laser intensity and electrostatic field strengths, the angle between them is continuously changed from to tt/2. The 
conclusions and outlook are provided in Sec. IIVI 



II. HAMILTONIAN OF AN ASYMMETRIC TOP MOLECULE IN THE PRESENCE OF THE FIELDS 

We consider a polar and polarizable asymmetric top molecule exposed to an homogeneous static electric field and 
a nonresonant linearly polarized laser. Our study is restricted to the regime of field strengths that significantly 
affects the rotational dynamics of the molecule, whereas its impact on the electronic and vibrational structure can be 
described by first order perturbation theory. We work within the Born-Oppenhcimer approximation, assuming that 
the rotational and vibrational dynamics can be adiabatically separated, and apply a rigid rotor description of the 
molecular systems. Furthermore, we neglect relativistic, fine and hyperfine interactions as well as couplings of different 
electronic states. In the laboratory fixed frame (LFF) {X,Y,Z), the Z-axis is chosen parallel to the polarization of 
the laser, and the direction of the homogenous electric field is taken forming an angle (3 with this axis and contained 
in the XZ plane. The molecular or body fixed frame (MFF) (x, y, z) is defined so that the permanent electric dipole 
moment is parallel to the z-axis, and for the considered systems the smallest moment of inertia is parallel to the 
X-axis. The relation between both frames is given by the Euler angles 57 = ((/>, 9, x) j43| . which are shown together 
with the field configurations in Fig. [TJ We only analyze molecules having the electric dipole moment parallel to one 
of the axis, and a diagonal polarizability tensor. Thus, the rigid rotor Hamiltonian reads 

H = Hr+Hs + HL, (1) 

where is the field-free Hamiltonian, and Hs and Hl stand for the interactions with the static and the laser field, 
respectively. 

In the absence of the fields, the rigid rotor Hamiltonian is given by 

Hr = BxJx + ByJy + BzJI (2) 
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where the angular momentum operators refer to the MFF, with Bi = ffi /2Iii being the rotational constant and In 
the moment of inertia around the principal axis of inertia i, with i = x,y and z. An estimation of the degree of 
asymmetry is provided by Ray's parameter, k ~ {2By — ~ Bx)/{By — B^), where the rotational constants of the 
considered system satisfy B^ > By > B^- For a symmetric top (Bx = By)i Ray's parameter takes the extreme values 
K = 1 and —1 for the oblate and prolate cases, respectively. 

For the regime of electrostatic field strengths Es considered here, we can neglect the interaction via the molecular 
polarizability, arriving at the following the Stark Hamiltonian 

Hs = -'Es-n=-Esficos9s, (3) 

where 9s is the angle between the permanent molecular electric dipole moment, fj, = i^iz, and the static electric field 
with cos 6*5 = cos /3 cos + sin /3 sin 6' cos and < /3 < 7r/2. 

Here, we consider a nonresonant laser field linearly polarized along the Z axis, EL(i) = i?niax5(i) cos(27rt^i)Z, with 
the frequency v, the field strength ii^max, and g{t) being the pulse envelope. We assume that is much shorter 
than the pulse duration or the rotational period, so that we can average over the rapid oscillations, which causes the 
coupling of this field with the permanent dipole moment to vanish [4J, |45|. In addition, we assume that the pulse 
duration is much longer than the rotational period of the molecular system, such that the states adiabatically follow 
the change of the field, and we restrict our analysis to the adiabatic limit g{t) — > 1. Hence, the interaction of the laser 
with the polarizability is the leading order term, and can be written as 

Hl = ^(a"^ cos^ e + ay"" sin^ 9 sin^ x) (4) 

2ceo 

with a-** = ajj — an, and an being the z-th diagonal element of the polarizability tensor, with i ^ x, y, z [46j . In this 
expression we have used (-Ef,) = I /ctQ, c being the speed of light, eo the dielectric constant, and expectation value 
indicates the time average. 

Our aim is to investigate the rotational spectrum of an asymmetric top molecule exposed to different field configu- 
rations. To do so, we solve the time- independent Schrodinger equation associated to the Hamiltonian ([T]), which even 
in the field-free case can not be solved analytically. Let us start by analyzing how the symmetries of this Hamiltonian 
change when the angle between the fields varies. 



A. Symmetries 

A detailed analysis of the symmetries of an asymmetric top rotor in the field-free case and exposed to a static 
electric field has been performed in Ref. [3^ . Here, we extend this study to the field configurations investigated in this 
work. 

The symmetries of the field-free Hamiltonian ([2]) are the spatial S0(3) rotation group and a subgroup of this group 
relevant to the symmetries in the presence of the fields is the molecular point group D2 = {E, C2,C2,C2}, where E 
stands for the identity and represents a rotation of tt around the i-axis of the MFF, with i = x, y, z. The action 
of these operators on the Euler angles are summarized in Table ID The -D2-group has four irreducible representations. 




FIG. 1: Laboratory and molecular fixed coordinate frames and field configuration. 
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TABLE I: Action of the symmetry operations on the Euler angles. 'The reflection can not be represented only by a rotation 
and the operation y — i- —y should be performed as well. 

The Wang states, defined as 

{JKMs}"" = -^{\JKM) + {-iy\J- KM)) , K>0 
I JO A/0)"' = |JOM), K = 0, 

with s = and 1, form the basis of these irreducible representations, characterized by the parity of J + s and K. 
The action of the elements of the L'2-group is C^lJKMs)"^ = (-1)^' | J/sTMs)"', with i = x,y,z and = J + K + s, 
Xy = J + s, and Xz = K. The states \ JKM) are the eigenfunctions of the field- free symmetric top rotor 

\JKM) = [-l)M-'<^^J-^Di,,^_An), (5) 

with Df[ k{^) being the Wigner matrix elements [43j . J the total angular momentum, and K and M the projections 
of J on the MFF z-axis and on the LFF Z-axis, respectively. To be self-contained, the definition and main properties 
of the Wigner matrix elements are in Appendix El 

For a field-free asymmetric top rotor, J and M are good quantum numbers, whereas, in contrast to a symmetric 
rotor, K is not well defined. For each M-value, there are four irreducible representations that depend on the parity 
of J -I- s and K . The eigenstates are degenerate with respect to M, and for a certain M and J, the corresponding 
eigenfunctions are linear combinations of Wang states \ JKMs)'^ with different K values. 

Since an external field defines a preferred direction in space, the symmetries of the corresponding Hamiltonian are 
reduced compared to the field-free case. As a consequence, the total angular momentum J is not a good quantum 
number, and only for certain field configurations M remains as a good quantum number. 

In the presence of a nonresonant laser field linearly polarized along the Z-axis, the symmetry operations of the 
Hamiltonian are the D2 point group, a rotation of an arbitrary angle 5 around the Z-axis Cz{5), a rotation of tt 
around an axis perpendicular to the Z-axis tilted at an angle a with respect to the X-axis C^^l"")! 8'iid the reflection 
in any plane including the Z-axis (this reflection is equivalent to flrst applying a twofold rotation around any axis in 
the XF-plane followed by the action of the operator C| or C| [131 )• Only M remains a good quantum number, and 
the levels with ±M are degenerate. Thus, we have eight irreducible representations for each \M\ > 0, but due to the 
reflection at the plane there is a twofold degeneracy and they can be effectively reduced to four representations being 
characterized by the parity of K and J -|- s, as in the field-free Hamiltonian. For M = 0, there are eight irreducible 
representations labeled by the parities of J, K and s. 

When an asymmetric top rotor is exposed to a static electric held (parallel to the Z-axis) or to both fields in the 
parallel configuration, i.e., /? = 0, the symmetry operations are C| from the D2 point group, the rotation Cz{5), and 
the reflection in any plane including the Z-axis. In these two cases, M is still a good quantum number, and the states 
M and —M are degenerate. For a certain |M|, the group has 4 irreducible representations characterized by the parity 
of K and the parity of s. Those representations with the same parity of K and M > are degenerate in energy. The 
symmetric top eigenfunctions ([5]), with defined parity of K , form a basis of this irreducible representation. For the 
M = case, these four representations are not energetically degenerate. 

For non parallel fields, M ceases to be a good quantum number. If the two fields are perpendicular, i.e., the electric 
fleld is parallel to the X-axis and /3 = 7r/2, the Hamiltonian commutes with only three symmetry operations: C|, a 
rotation of tt around the X-axis (tt), and the reflection axz on the XZ-plane where the fields are contained, as 
well as their combinations. Using the field- free symmetric rotor wave functions, we can construct a basis for crxz, 



iJKMq)" = —={\JKM) + (-1)«|J- A' - il/)) 
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with M and/or if 7^ 0, and 

I JOOO)'" = I JOO) 

with M = K = 0, where g = and 1, and (Jxz\JKMqY = {-lY^+^+'i\J KMq)" ; and one for Cx(7r) 



with M ^ 0, and 



\JKOO)^ ^ \JKO), 



with M — 0, wherep = and 1, and satisfying Cx{'!^)\J KM-p)^ = {—1)'-'+p\JK Mp)-^ . The group of all the symmetry 
operators, C|, axz and Cx(7r), has eight different irreducible representations according to the parity of M + K + q, 
J + p and K. A basis of these irreducible representations is 

\JKMqp)^ = ]^{\JKM) + {-iy\J -K-M) 



with M 7^ and if 7^ 0, 



with A/ 7^ and if = 



+ {-lY\JK - M) + J - KM)), 



I JOMOp)^ = ^(I^OM) + JO - M)), 



I JA'OqO) . = -^i\Jm + (-l)'l^ - ^'0»: 



with A/ = and if 7^ and 

IJOOOO)- = I JOO), 

with K = A/ = 0, where q^Q and 1, p = and 1, and the parity oi M -\- K + q, K and J + p are preserved. 

Finally, when the fields form an angle < /? < 7r/2, the Hamiltonian is invariant under the reflection axz and the 
rotation C|. We have, therefore, four irreducible representations depending on the parity of Af + K + q and of 
and the corresponding basis is {\JKMq)'^}. 

For an asymmetric top rotor exposed to any of these field configurations, the field-dressed spectrum exhibits many 
avoided crossings between energetically adjacent states of the same symmetry. When the spectrum is analyzed as the 
strength of one of these fields is varied or the angle between them, these avoided crossings should be distinguished 
from the real crossings taking place between levels of different symmetry. Hence, we solve the Schrodinger equation, 
by expanding the rotational wave function in a basis that respects the symmetries of the corresponding Hamiltonian. 
As a consequence, the coefficients of these expansions fulfill the properties of the basis vectors of the corresponding 
irreducible representation. For computational reasons, we have cut the (in principle) infinite scrie to a finite one 
including only those functions with J < Jmax, and for a certain J all (2 J + l)-values of K, and, analogously, for M 
in the case < /3 < 7r/2. The size of the Hamiltonian matrix increases as J^^ax and J^ax fo^' the < /3 < 7r/2 and 
P = configurations, respectively. In this study, we have used Jmax = 24, and the convergence is reached for the 
states analyzed here. Several matrix elements are presented in Appendix 1X1 

The field- free states are labeled by the notation Jk^.K-M, where Ka and Kc arc the values of K on the limiting 



symmetric top rotor prolate and oblate cases, respectively |47| . For reasons of addressability, we use this notation 
for the field-dressed states, even if J and/or M are not good quantum numbers. Thus, JKaKc_M refers to the level 
that is adiabatically connected as /, Es and/or (3 are modified with the field-free state Jk^k^M. The irreducible 
representation to which the states belong is also indicated. Analogously to a symmetric top molecule exposed to 
combined fields [IHl , the final labels of the states depend on the path followed on the parameters to reach a certain 
field configuration, i.e., monodromy is observed. Since each interaction breaks different symmetries of the field-free 
Hamiltonian, the order the fields are turned on determines the evolution of the field-dressed states. The complexity 
of the spectrum is characterized by the amount of genuine and avoided crossings among the states, the symmetry of 
the two levels determine the type of crossing that they may suffer as one of the field parameters [Es, I or /3) is varied, 
and, therefore, if the corresponding labels may or not be interchanged. 



6 



III. RESULTS 

In this section, we illustrate the impact of the external fields on two asymmetric top molecules: fluorobenzene 
(CgHsF) and pyridazine (C4H4N2). Their data are summarized in Table IIIl according to Rcfs. [isl - fs^ and their 
structure is shown in Fig. [2] They are characterized by a different degree of asymmetry, the fluorobenzene is 
intermediate-prolate with k = —0.5879, and the pyridazine is near-oblate with k = 0.8824. The permanent dipole 
moment of pyridazine is around 2.5 times larger than in fluorobenzene. The asymmetry of the polarizability tensor 
is very similar for both systems, for fluorobenzene, a'-^ = 4.298A^ and a^^ = 3.848A^, whereas for pyridazine 
^zx 4 52A and a^^ = 4.45A . Since the rotational constants of pyridazine are larger than for fluorobenzene, for 
the same laser intensity, a weaker impact on the former should be expected. In the following, we carry out a study 
of the spectrum of these two systems as the parameters that characterize the fleld configurations, Es, I or /3, are 
modified. For the sake of simplicity, we analyze the energy, the orientation, the alignment and hybridization of the 
states which belong to the representation with K and M + K + q even (if /3 7^ 7r/2), and J + p even (if f3 = tt/2). 
Note that they represent well the main physical features observed in the overall spectrum, and similar behavior and 
properties are therefore obtained for the other representations. 




FIG. 2: Structure of the (a) fluorobenzene and (b) pyridazine molecules. 



Fluorobenzene Pyridazine 



(MHz) 


1716.916 


3055.485 


By (MHz) 


2570.624 


6048.613 


B^ (MHz) 


5663.72 


6235.680 


K 


-0.5879 


0.8824 


ti. (D) 


1.66 


4.14 


Ctxx (A ) 


7.141 


5.84 


ayy (A ) 


10.89 


10.29 


a.z (A') 


11.439 


10.35 



TABLE II: Relevant data for fluorobenzene [i^ - lsot and pyridazine [5ll.[53|. 



A. Impact of a linearly polarized laser field 

The dynamics of the molecule in the presence of a linearly polarized laser field depends strongly on the anisotropy 
of the polarizability tensor. The interaction with the laser ^ exhibits several critical points. The minimum value 
Hi^ = —Ia^^/2eoc is reached for 9 = or n and any value of x- The interaction achieves three maxima at 9 = tt/2 and 
X = 0, TT or 27r satisfying that Hl — 0, and two saddle points at 6* = 7r/2 and x = '^/'^ or 37r/2 with = — /a''^/2eoc. 
All of them are shown in Figs, ^a) and (b), where Hl is plotted for x = '^/'^ and < 9 < 27r, and for 9 = Tr/2 
and < X < Stt, respectively. The main difference between the two systems is that the values for Hl at the saddle 
point 9 = Tr/2 and x = ''^/'^ is smaller for pyridazine compared to fluorobenzene, and the fact that the shape of 
Hl as a function of 9 and for fixed x is significantly flatter for the former. Taking into account that the rotational 
constants arc much larger for pyradizine compared to fluorobenzene we observe that for the same laser intensity /, the 
pyridazine wave function is more widespread with respect to 9 and x than for fiuorobenzene. While the minima are 
responsible for the molecular alignment, the maxima or saddle points correspond to an " antialigned" wave function, 
that is the dipole moment points perpendicular to the field direction. 
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FIG. 4: Ground state energy (a), (cos'^ 0) (b), (sin^ x) (c), and {K'^) (d), as a function of the intensity of a linearly polarized 
laser field, for the fluorobenzene (solid) and pyridazine (dash) molecules. 

To get a better physical insight into this interaction, we present its impact on the ground state energy, and the 
expectation values (cos^ O), (siv? x)^ and (K"^) in Figs.|4|a), (b), (c), and (d), respectively. While, for both molecules, 
the energy as a function of / shows a similar decreasing behavior, their values being indistinguishable on the scale 
of Fig. m^a), significant differences are observed for the other quantities. In the very strong laser field regime, the 
probability density of the fluorobenzene ground state tends to concentrate around the minima, and one should expect 
that (cos^fl) — > 1, and (sin^ x) ~^ 0-5 for very large intensities. This last relation holds because the x-coordinate 
does not play any role in the absolute minima of Hl, and one concludes that the probability density should be 
uniformly distributed with respect to x- Numerically we obtain that for the fluorobenzene ground state, (cos^ (?) 
increases until 0.76 for / = 10^^ Wcm~^, decreasing thereafter, and (sin^ x) reaches a plateau with a constant value 
0.58 for / > 4.2 • 10^" Wcm~^. The barrier height of Hl as a function of 6 for a certain value of x, see Fig. ^a), is 
around 9.2 times smaller for pyridazine than for fluorobenzene, and for the former, the rotational constants are larger, 
whereas the polarizability anisotropics a^^ and a^^ are of the same order for both molecules. Hence, compared to 
fluorobenzene, the pyridazine ground state wave function should be spatially stronger delocalize with respect to 9, 
and, therefore, less aligned for the same laser intensity. We obtain here ^cos^6') ~ 0.49 for / = 10^^ Wcm^'^. As a 
consequence, the spreading of the wave function for x is not observed in pyridazine: ^sin^ x) shows a maximum of 
0.75 for / = 3 • 10^^ Wcm~^ and slightly decreases with further increasing /. For the field-free ground state, we have 
(K^) = 0, and as / is increased (K^) follows a similar evolution as (sin^x)- Foi' fluorobenzene, (K^) achieves the 
value 0.14 for / « 7.2 x 10^° Wcm~^, followed by a plateau-like behavior around {K'^} ~ 0.18 for larger intensities 
for fluorobenzene, while it keeps an increasing trend up to 1.38 for / = 5.9 x 10^^ Wcm~^ and decreasing smoothly 
afterwards for pyridazine. Let us emphasize that, for a field-free asymmetric rotor, K is not a good quantum number, 
and an eigenstate already shows a certain amount of JiT-mixing, but in the strong laser field regime, the second term 
of the laser interaction Q should impact and enhance this i^-mixing. 
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FIG. 5: (a)-(c) Energies and the expectation values (d)-(f) (cos^s), (g)-(i) (cos^ 6*) and (j)-(l) (cos^ x) for a constant field 
Es = 20 kVcm^^ as a function of the intensity of the laser field for /3 = 7r/6, 7r/4 and n/3 for the first states with both 
M + q + K and K even for fluorobenzene. The states are OooO (solid black), loil (solid green), loiO (dash dotted black), 2o22 
(dash dotted green), 2o2l (dash black) and 3o32 (dash green). The spectrum (a, b, c) contains also highly excited states (very 
thin lines). 

B. Constant static electric field and increasing laser intensity 

In the presence of an additional static field, the interaction is given by Hs + H^, see eqs. ([3]) and ([4]), and the 
dynamics is significantly more complicated. The amount of extremal points of this potential and their character 
strongly depend on the field parameters as well as on the molecular polarizability and permanent dipole moment. 

For fluorobenzene and pyridazine, we show in Figs. [5] and [6] the dependence of the energies (panels (a), (b), and 
(c)), expectation values {cos9s) (panels (d), (e), and (f)), (cos^ 6*) (panels (g), (h), and (i)), and (cos^ x) (panels (j), 
(k), and (1)), as functions of the laser intensity for a constant electric field Es — 20 kVcm^^, and /? = 7r/6, 7r/4 and 
7r/3, respectively. For the sake of simplicity and without loosing generality, we restrict our analysis to the energetically 
lowest-lying six states for the irreducible representations with both K and M + K + q being even, which are the levels 
Jk^,k^M = OooO, loiO, loil, 2o2l, 2o22 and 3o32 for fiuorobenzene, and Jk,,k^M = OooO, loiO, loil, 2o2l, 2o22, 
and 22i2, for pyridazine. To illustrate the complexity of the spectrum, we have included in the energy panels highly 
excited levels with the same symmetry (very thin lines). The adiabatic following has been done by increasing first 
the strength of the static electric field up to 20 kVem~^, the static field being tilted by an angle /3 with respect to the 
Z-axis. This is done for / = and yields the labeling of the states in the presence of the static field. Thereafter, the 
laser intensity is increased. For both molecular systems, these levels are high-field-seekers, and their energies decrease 
as / is increased. For a given laser intensity, the lowering in energy (compared to the field-free value) increases with 
decreasing angle between both fields. Since all the states included in these figures possess the same symmetry, we 
encounter exclusively avoided crossings of energetically adjacent states. We have assumed that the avoided crossings 
are traversed adiabatically as / is increased (according to Landau-Zener transition theory), and consequently the 
character of the involved states is interchanged. These avoided crossings, which are not distinguishable on the energy 
scale panels (a)-(c) in Figs. [S] and [51 strongly affect the orientation and alignment features of these levels. 

Let us start analyzing the results for fiuorobenzene in Fig. [SJ The orientation of the corresponding wave functions 
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is illustrated by the expectation value (cos 0s), with 0$ being the angle between the static electric field and the 
molecular fixed z-axis, which coincides with the direction of the permanent dipole moment. Only the states OpoO and 
3o32 present a significant orientation with respect to the static electric field direction, that is reduced as /3 is increased. 
The ground state satisfies {cos 9s) > 0.70 for the three /? values, and it has a plateau like behavior, with a minor 
positive or negative slope as / is increased. The numerous avoided crossings have significant impact on the other 
levels especially for /3 = tt/6 and 7r/4, and the underlying states might evolve from a strongly oriented configuration 
into a weakly oriented or antioricnted one. As an example, the state 3o32 after suffering for / w 1.7 • 10^° Wcm~^ an 
avoided crossing with the non-oriented 2o2l state looses its strong orientation. By further increasing / the 2o2l level 
suffers another avoided crossing, which provokes a local maximum in {cos 6s) for / w 2.4 • 10^° Wcm~^. Since for 
/3 = 7r/3 the width of all the avoided crossings is larger for all the levels, the orientation shows a smooth evolution as 
/ is varied. A nonresonant linearly polarized laser field provokes the alignment of the wave function along the Z-axis, 
i.e., (cos^ 0) tends to increase and ultimately approach the value 1 as / is increased (Fig. [5](g, h, i)). This process 
competes now with the orientation due to the static field. Only the ground state alignment keeps an increasing trend 
as / is enhanced, and for the given three configurations (cos^ 0) > 0.75 for / = 10^^ Wcm^^. An additional electric 
field at an angle /3 = 7r/6 or tt/4 favors the alignment and ^cos^ 6*) is larger than without static field, sec Fig. Hljb), 
whereas as the angle between the fields is augmented the values achieved for (cos^ 9) come closer to those of Fig. [41(b). 
The level Igil does not achieve a large alignment, and for (3 = 7r/4, (cos^ 6*) exhibits a broad well that for f3 = tt/S 
is even wider since the coupling between the states changes. Around the avoided crossings, the wave function of the 
involved states alternate regions of significant alignment with other characterized by broad distribution as the laser 
intensity is varied, e. g. see the aligment of the states 3o32 and 2o2l for / « 1.7 • 10^° and 3.5 • 10^° Wcm~^ in Fig. [5] 
(g) and (h); respectively. The impact of the avoided crossings is also noticeable for f3 = n/3. The behavior of (cos^ x) 
strongly depends on the considered level. Since the molecules exhibit a strong orientation with the electric field, the 
contribution in x should be to increase the term (sin'^ x) which will give rise to a decrease of the energy. For the 




2.5 5 7.5 10 2.5 5 7.5 10 2.5 5 7.5 10 

Wcm2] /[lOi" Wcm2] /[iQio Wcm^] 



FIG. 6: (a)-(c) Energies and the expectation values (d)-(f) (cos^s), (g)-(i) (cos^ 6*) and (j)-(l) (cos^ x) for a constant field 
Es = 20 kVcm^^ as a function of the intensity of the laser field for /3 = 7r/6, 7r/4 and 7r/3 for the first states with both 
M + q + K and K even for pyridazine. The states are OooO (solid black), loil (solid green), 2o22 (dash dotted black), loiO 
(dash dotted green), 2o2l (dash black) and 22i2 (dash green). The spectrum (a, b, c) contains also highly excited states (very 
thin lines). 
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states OqoO and loil, (cos^ x) decreases as / is increased, and for both states there exists some region where (cos^ x) 
keeps a smooth behavior. For the other states, this expectation value is also affected by the presence of avoided 
crossings, and (cos^ x) alternates between increasing and decreasing behavior as a function of the laser intensity. 
Since the pyridazine possesses a larger permanent dipole moment than fluorobenzene, the impact of the static field 
is larger, and also dominates the dynamics, see Fig. [6] For the ground state and first excited one, (cos 9s) is only 
weakly affected by the laser field, for the former {cos 9s) ~ 0.75 independently of / and /3, while for the loil level 
{cos9s) > 0.5. The other levels present a right-way orientation, that can be converted from a strong to a mild one 
when an avoided crossing is encountered. Regarding the alignment, for most of the states ^cos^ 9'^ shows a smooth 
behavior as / is enhanced for the three values of /3. The impact of the avoided crossings on this expectation value is 
not very pronounced because most of the states show a weak alignment with a similar value of (cos^ 0). Analogously 
to fluorobenzene, the ground state has a larger alignment for /3 = 7r/6 and 7r/4 than in the absence of the static field, 
while for /? = 7r/3 other states are found with larger alignment. For stronger fields, (cos^ x) monotonically decreases 
as / is enhanced for all the states, and the slope is more pronounced compared to fiuorobenzene. 



C. Constant laser intensity and increasing electric field strength 



The pendular limit of an asymmetric top molecule in the presence of a strong electrostatic field was investigated 
by Kanya and Ohshima [U using a power series expansion in fiEs- Their analytical expression for the energy, which 
neglects the contribution of terms in powers equal or smaller than {ijEs)~^^^ (see eq. (26) in Ref. Isil) allowed us for a 
straightforward comparison to our numerical calculations. For Es = 100 kVcm~^, the energy and orientation cosine 
of the ground state agree within 0.05% and 0.025% for fiuorobenzene and 0.04% and 0.019% pyridazine, respectively. 
Note that for highly excited states these relative errors increase. 

For a constant laser field / = 10^" Wcm~^, we now investigate the impact of increasing static field strength for three 
different configurations. Again, we consider the energetically lowest-lying six states with the irreducible representation 
for K and M + K + q being even. These levels have been adiabatically followed as the laser intensity is raised from 
/ = to 10^° Wcm"2 for Eg = 0, we label them, and finally the electrostatic field is turned on forming an angle 
/3 and its strength is increased. Thus, for both molecules the levels are OqoO, IqiO, IqiI, 2o20, 2o2l, and 2o22. For 
fiuorobenzene and pyridazine, we present in Figs. [7]and[8l the evolution of the energies (panels (a), (b), and (c)), 
(cos^s) (panels (d), (e), and (f)), and (cos^ 9) (panels (g), (h), and (i)), as Es is increased, for / = 10^" Wcm~^, and 
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FIG. 7: (a)-(c) Energies and the expectation values (d)-(f) {cos 9s), (g)-(i) (cos'^ 6) for a constant I = 10^'^ Wcm^'^ as a function 
of the strength of the static field for /3 = tt/6, 7r/4 and tt/3 for the energetically lowest states with both M + q + K and K even 
for fluorobenzene. The states are OooO (solid black), loiO (solid green), loil (dash dotted black), 2o20 (dash dotted green), 2o2l 
(dash black) and 2o22 (dash green). The spectrum (a, b, c) contains also highly excited states (very thin lines). 
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/3 = 7r/6,7r/4 and 7r/3, respectively. We remark that for such a weak laser intensity / = 10"'^'^ Wcm~^ a static electric 
field of Eg > 17.08 and 7.26 kVcrn"^ for fluorobenzene and pyridazine, respectively, provides the larger contribution 
to the external field Hamiltonian Hs + Hl . If the interaction with the static field is dominant and much larger than 
the laser one, i.e., Egfi ^ la^^ /2eoc with i = y or z, the absolute minima of the potential Hs + Hl are at (0, /3,x), 
with ^ 0, 27r and x ~^ ""/S, 37r/2. Thus, the wave function will be oriented toward the electric field direction with 
(cos^ 0) — cos^ /? and {cos 8s) 1. However, these extremal points do not provoke any effect on the x coordinate, 
because in the strong static field regime, the variation of Hs + H^ as a function of x represents a very shallow 
minimum and does not give rise to a localization of the wave function. 

In the weak field regime there are, for both molecules, several levels in addition to the ground state that are high- 
field seekers, and remaining levels are for low values of Es low- field seekers. These levels present a mild wrong- way 
orientation with (cos^s) < 0. Since in this regime the slope of the variation of the energy with Es might be positive 
or negative, this favors the presence of sharp avoided crossings. As a consequence, the orientation of a pair of levels 
involved in an avoided crossings, i.e., (cos 6*5), suffers drastic variations over tiny ranges of the field strength. As the 
field strength is increased, we encounter the pendular regime: all the states are high- field-seekers, and they are strongly 
oriented along the static field direction. Indeed, for Es = 100 kVcm"-'^, we have that (cos 6*5) > 0.60 and 0.70 for the 
considered fluorobenzene and pyridazine levels, respectively. In this regime, we still encounter avoided crossings but 
they arc much wider. Due to the competition between both fields, these states do not achieve a significant alignment, 
see panels (g)-(i) in Fig. [7]and|Sl We observe that (cos^ 6') approaches cos^ /? in the strong field regime. 



D. Orientation and 2-D alignment by means of perpendicular fields 



External fields provide a tool to control the molecular dynamics: specifically it has been shown that an elliptically 
polarized laser allows for 3-D alignment of asymmetric molecules, i.e., the system is aligned in all spatial directions [2y, 
El, m, [s^l ■ We show here that the combination of an electrostatic field with the linearly polarized laser gives rise to 
orientation in one direction and alignment in the other two. We hereby focus on the case /? = 7r/2. In the strong 
static field regime, we have 9 — 7r/2 and ^ — >■ 0, whereas, in order to decrease the energy the term sin^ x in the laser 
Hamiltonian ^ should increase and approach one, i.e., x ~^ ""/S, 37r/2. Then, the molecule will be fixed in space, and 
the wave function should be concentrated in the proximity oi 9 ^ ''■/2, x ~ '''/2, 37r/2 and ~ 0, 27r. To demonstrate 
this behavior we have computed the ID probability density distribution in each one of the three Euler angles by 
integrating the square of the wave function in the other two angles. In Fig. E^a), (b), and (c), these probability 
density distributions 0(0) sin 0, $(</>) and S(x) are plotted as a function of 0, and X: respectively, for the ground 
state of pyridazine interacting with the orthogonal fields for different strengths. The probability density distributions 




FIG. 8: Same as Fig. [7| but for pyridazine. 
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FIG. 9: One dimensional probability density distribution in each Euler angle for the ground state of pyridazine for Es = 20 
kVcm"^ and / = lO" Wcm^^ (solid), Es = 50 kVcm^^ and / = lO" Wcm^^ (dash dotted), Es = 20 kVcm^^ and I = 5 - lO" 
Wcm"^ (dash). 

for Es = 20 kVcm^-'^ and / = 10^^ Wcm~^ show the expected behavior, and this state presents a significant orientation 
along the X-axis of the LFF with {cos9s) = 0.749 and (cos^ 0) = 0.191, whereas the alignment in the other two 
Euler angles is also pronounced and we get (cos^ x) = 0.143 and (cos^ (/>) = 0.661. A larger orientation is achieved 
if the static field strength is enhanced to Es = 50 kVcm~^ keeping the same laser intensity: (cosOs) = 0.849 and 
(cos^ 0) = 0.118, for the distribution in x and and we find (cos^ x) = 0.122 and (cos(/>) = 0.723, respectively. An 
opposite effect has the enhancement of the laser intensity to / = 5 • 10^"'^ Wcm~^ and keeping Es = 20 kVcm~^. The 
orientation in 9 is slightly reduced {cos 9s) = 0.729 and (cos^ 0) = 0.222, but (cos^ 0) increases only to 0.682, that is, 
the laser does not affect the azimuthal angle and in Fig. |9|b) the distribution $(0) is indistinguishable compared to 
the corresponding one for Es = 20 kVcm~^ and / = 10^^ Wcm~^. The alignment with respect to x becomes stronger 
(cos^ x) = 0.075. Of course, one should keep in mind that this effect is very sensitive to the fields strengths and their 
configuration, as well as and the molecular parameters, i. e., the permanent dipole moment, polarizability anisotropics 
and rotational constants. For fluorobenzene, a similar phenomenon could be found, but the wave function shows a 
slightly less pronounced orientation and alignment. If the fields are not perpendicular, the competition between both 
interactions will reduce the alignment achieved in the (j) a-nd x angles. 

E. Influence of the inclination of the fields 

The inclination of the fields plays an essential role on their impact on the rotational dynamics of the system, and 
as already discussed, the symmetries are drastically modified as /3 is varied. Indeed, /? is another parameter in the 
Hamiltonian that enriches the physical phenomena observed, and its variation can provoke the appearance of avoided 
crossings between energetically adjacent states of the same symmetry. For pyridazine, we represent in Figs. [10] (a) 
the energy, (b) {coa9s), (c) (cos^0), and (d) (M"^), as a function of /3 with Es = 20 kVcm~^ and / = 10^^ Wcm^^, 
for the previous set of states with K and M + K + q even and following the same labeling as in Sec. IIII Bl For 
these field strengths, the static electric field interaction dominates over the laser field interaction, similar to Fig. |51 In 
general, the energies show a smooth behavior as /? is varied, and depending on the state they exhibit an increasing or 
decreasing trend, see Fig. lTUf a). e.g., the ground state energy increases from —5 cm~^ to —4.6 cm"-'^ from /3 = to tt/2, 
respectively. The energy gap between the first four states is large enough to prevent the presence of avoided crossings 
among them, the first avoided crossing being between the fourth and fifth excited states, 2o2l and 22i2, for f3 « 1.096 
(close to 37r/8). Regarding the orientation and alignment along the electric and laser fields directions, respectively, 
different behaviors are observed. The ground state keeps a significant and approximately constant orientation with 
(cos 0s) > 0.75 for any value of /3, as the electric field is rotated away from the Z-axis the ground state probability 
density follows this field. In contrast to this, since the laser interaction is not dominant, its alignment is drastically 
reduced from (cos^ 0) ~ 0.75 to 0.22 when /3 increases from to 7r/2. For any field configuration, the Igil level shows 
a moderate orientation and alignment with a plateau-like behavior for {cos 9s) and (cos^ 0). Compared to the ground 
state, the 2o22 level presents a similar orientation and alignment for /3 = 0, but a very different evolution of these 
features as /3 is varied, and for /? = 7r/2, it keeps a moderate alignment and a weak orientation. For parallel fields, M 
is a good quantum number, and a non-parallel configuration allows the interaction and mixing between states with 
different field-free M-value. This phenomenon is illustrated by means of the expectation value {M'^) in Fig. [TCT d). 
As the angle /3 is increased, the evolution of ^M^) strongly depends on the character of the corresponding level. For 
the ground state, (M^) increases as f3 is enhanced, and for (3 = tt/2, it reads (M^) = 0.746. The (M^) value of the 
loil level is close to 1 for /3 < tt/8, but for larger values of /3 decreases to 0.6 for /3 = tt/2. In contrast, for the other 
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FIG. 10: (a) Energies and expectation values (b) (cos6g), (c) {cos^ 9) and (d) (A/^) for pyridazine in the presence of a static 
field Es = 20 kVcm~^ and a laser field I = 10^^ Wcm"'^ as a function of /? for the energetically lowest states with both 
M + q + K an K even. The states are OooO (solid black), foil (solid green), 2o22 (dash dotted black), loiO (dash dotted green), 
2o2l (dash black) and 22i2 (dash green). 




FIG. 11: Width AE (a) and electrostatic field strength Es (b) at the avoided crossing taking place between the states 2o20 
and 2o2l for pyridazine, for / = 10^° Wcm~^ and different inclination angles /3. 

analyzed level vifith M = 1, 2o2l, the interaction with states with larger M is dominant for (3 > 7r/4, and (M^) ~ 3 
for P = 7r/2. For the considered M = 2 levels, 2o22 and 22i2, the mixing with states with lower M is dominant, and 
(M^) is smaller than in the parallel configuration, e.g., for the 2o22 level and (3 = 7r/2, we have (M2) = 0.4063. 

With varying inclination angle /3 between both fields the avoided crossings leave their fingerprints in the relevant 
observable. For pyridazine, the states 2o20 and 2o2l belong to different irreducible representations for /3 = and 
7r/2, and to the same one for < /3 < tt/2. (Note that the labeling of the states has been done in the same way 
as in Fig. [5]). For non collinear fields, they suffer an avoided crossing which we have traced for / = 10^*^ Wcm^^ 
and different values of /3 in Fig. [TT] The results for the minimal energetical width AE = |£'2q20 — ^2o2i| the 
electrostatic field strength at which this minimum appears are presented in Figs. II If a) and (b), respectively, as the 
angle /3 is varied. For perpendicular fields, the levels suffer a real crossing and are accidently degenerate AE ~ for 
I = 10^° Wcm~^ and Es = 39.78 kVcm~^, whereas for /3 = 0, they possess a different magnetic quantum number M 
and exhibit a symmetry-related crossing but now for Es = 47.55 kVcm~^. As we see in panel (a), AE increases till 
it reaches the maximal value of 5.33 • 10~^ cm~^ for (3 ~ 7r/4. decreasing afterwards to for /? = 7r/2. The static field 
strength at which the avoided crossing takes places, see panel (b), decreases monotonously as (3 increases. Indeed, 
the variation of Eg with (3 is well-matched by the following function 3.88 cos(2.06/3) -I- 43.67, and it is reduced by 7.74 
kVcm~^ when (3 goes from to tt/2. 

IV. CONCLUSIONS 

In this work, we have investigated the impact of a combination of an electrostatic and a nonresonant linearly 
polarized laser field on the rotational spectrum of asymmetric top molecules. This study has been performed in 
the framework of the Born-Oppenheimer approximation considering that the vibrational and electronic dynamics are 
not affected by the external fields. Our analysis is restricted to a rigid rotor description of molecules having their 
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permanent dipole moment parallel to one axis of inertia, and the polarizability tensor is diagonal in the basis formed 
by the principal axis of inertia. We have analyzed the symmetries and irreducible representations of the Hamiltonian 
for the different field configurations. Numerically, each irreducible representation has been treated independently, 
by expanding the wave function in a basis respecting the corresponding symmetries. This procedure allows us to 
distinguish the avoided crossings from genuine ones involving states of the same and different symmetry, respectively. 
The presence of the avoided crossings in the field-dressed spectrum affects the directional properties of the molecule, 
they might significantly alter the spectroscopy as well as the sterodynamic of the system. The labeling of a certain 
state and the passage through the avoided crossing that it suffers depend on the way the symmetries are broken, 
i.e., on the temporal sequence followed to turn on the fields, which should be taken into account to determine the 
adiabaticity of a certain process. 

The richness and variety of the resulting field-dressed rotational dynamics have been illustrated by analyzing 
the energetic shifts, as well as the orientation, alignment and the hybridization of the angular motion. As prototype 
example we have investigated the fluorobenzene and pyridazine molecules. For three field configurations, the evolution 
of a certain set of states belonging to a certain symmetry has been analyzed with varying electrostatic field strength 
or laser intensity. Different types of behavior were observed, depending on the dominant field interaction as well as 
on the considered molecular system, through its rotational constant, dipole moment and polarizability tensor. Due 
to the competition between both interactions, the features of the rotational spectrum are significantly changed as 
the field parameters arc modified. In the strong laser field regime, the presence of an intense electric field reduces 
the orientation of the ground state, especially as /3 is rotated from zero to larger values, and highly excited levels 
only present a very weak alignment. Whereas, if the electrostatic field is dominant the states are oriented along 
its direction, and they only present a mild alignment along the Z-axis in the LFF. We have shown that a proper 
combination of non-collinear fields gives rise to a strong orientation along the static field direction together with a 
2D-alignment on the other two axis of the molecule. In particular for f3 = 7r/2, the molecular plane is fixed onto the 
Xy-plane of the LFF, and the orientation is along the X axis. Finally, we have also investigated the role played by 
the inclination angle of the fields /3, by analyzing the spectral properties of several states, the loss of the azimuthal 
symmetry has been quantified by the expectation value (M^), which is a conserved magnitude for parallel fields. An 
avoided crossing between two states has been traced as /3 is modified, the electrostatic field strength at which it takes 
places varies within a few kVcm^^, and the corresponding energetical width will allow us to compute the adiabaticity 
of the crossing once the variation of the field strength is known. A natural extension to this work would be to consider 
other molecular systems, especially different conformers of the same molecule, looking for specific phenomena that 
might help to distinguish between the molecules. 
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Appendix A: Wigner and Hamiltonian matrix elements 

The field-free eigenstates of a symmetric top molecule, see eq. (O, are proportional to the Wigner matrix elements, 
D'Ij j^{Q), which are defined as 

Di^i,{n) ^ e-^''Uij^^{0)e~^^^, (Al) 

where dfj ^ (6) are the reduced Wigner matrix elements [1^ . To evaluate the matrix elements of the Hamiltonian, we 
have used the following properties of the Wigner matrix, the complex conjugate 

DilA^) = {~ir-''Di,,^_^{n), (A2) 
the relation between the Wigner reduced matrix elements reads 

= i-^)"' <^,„Ae) (A3) 

<',n{0) = (-l)"'""di„.,-„(^) (A4) 
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and the integral of the triple product of Wigner matrices 



Sir' 



Ji J2 J \ Ji J2 J 
Ml M2 -Ml \Ki K2 -K 



(A6) 



where 



are the 3J Symbols. 



Jl J2 J3 

For completeness, we provide the non-zero matrix elements appearing in the evaluation of the full Hamiltonian. 
For the field- free Hamiltonian Hr we have 



{J'K'M'\Hr\JKM) = {AJ{J + 1) + CK^)5j,j6k',kS_ 



M'M 



+ By/J{j + 1) - K{K + l)y/J{J +1)- [K + 1)(A' + 2)5j, 

where A, B and C are defined as A = {B^ + By)/2, B = {B^ - By)/A, and C = (25^ - B^^ - By)/2 [H. The 
Stark interaction Hs ^ rewritten in terms of Wigner matrix elements, is given by 



(A7) 



Hs = -fiEscosOs = -^lEshs = -fiEg { cos pDl^m + sinp^/ -{Dl^^oi^) - DloiQ)) 



Thus, the non-zero matrix elements are 

(JKM 



(J KM + 1 



(J + IKM 



(J + IKM + 1 



hs\JKM) 
hs\JKM) 



MK 
J(J + 1) 

K sin /3 



2J{J + 1) 
hs\JKM) = cos/3 



cos/3 

^J{J+l)-M[M + l) 



[(J+1)2-M2][(J+1)2-X2] 



hs\JKM) = -'■ 



y (2J + 3)(2J+1)(J+1)2 
in/3v/[( J + 1)2 - iC2]( J + M + 1){J + M + 2) 



(J + IKM - 1 
The laser Hamiltonian Hl jS]) takes now the form 

a^x ^ ^yx 



2(J-Hl)v/(2J + 3)(2J+1) 
hs\JKM) = -{J + 1K - M + l\hs\JK - M) 



I I 

hL = - 



2eoc 



2enc 



g D^n) - '^[DUn) + Dl_^m + 



yx 



(A8) 



In this expression we have omitted the terms which only introduce a shift in the energy. Performing the integrals 
corresponding to each term, we get 



{JKM\hL\JKM) = 



^zxj^^zy\ [3Af2 - J( J -h l)][3/s:2 _ 1)] 



{J +lKM\hL\JKM) = (a^^ + a^2^) 



{J + 2KM\hL\JKM) = 



zx I 



2J(J-H1)(2J- l)(2J + 3) 

MKy/[{J + 1)2 - A/2] [(J + 1)2 - i^2] 

J( J + 1)( J -I- 2)^(2 J + 3)(2 J + 1) 

V[(J+2)2-if2][(J+i)2_^2] 



v/[(J+2)2-M2][(J+l)2-/\/2] 



(J + 1)(J + 2)(2 J + 3)V(2J+ l)(2J + 5) 
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{JKM\hL\JK + 2M) 



a^^[3M2- J(J+1)] 



X 



- (A- + l)2]( J _ K){J + K + 2) 
2J(J + 1)(2J + 3)(2J- 1) 



{J +lKM\hL\JK + 2M) 



X 



^[{J - K)^ - l]iJ - K){J + K + 2) 
2J(J + 1)( J + 2)^(2 J + 1)(2J + 3) 



(J + 2ii:A/|/ii|Jif + 2A/) 



-a«^\/[(J + 2)2 - M2][( J + 1)2 _ 1/2] 



v/[(J-A02-i](j-/v)(J-i^^ 



{j'K'M'\hL\JKM) 



4( J + 1)( J + 2)(2 J + 3)V(2J+l)(2J + 5) 
(J' - K' - M'\hL\J -K- M). 
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